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ABSTRACT 

We perform the largest currently available set of direct N-body calculations of young star 
cluster models to study the dynamical influence, especially through the ejections of the most 
massive star in the cluster, on the current relation between the maximum-stellar-mass and the 
star-cluster-mass. We vary several initial parameters such as the initial half-mass radius of the 
cluster, the initial binary fraction, and the degree of initial mass segregation. Two different 
pairing methods are used to construct massive binaries for more realistic initial conditions 
of massive binaries. We find that lower mass clusters 10^^ do not shoot out their 
heaviest star. In the case of massive clusters 1000 Mq), no most-massive star escapes 
the cluster within 3 Myr regardless of the initial conditions if clusters have initial half-mass 
radii, ro.5, ^ 0.8 pc. However, a few of the initially smaller sized clusters (ro.5 = 0.3 pc), 
which have a higher density, eject their most massive star within 3 Myr If clusters form with 
a compact size and their massive stars are born in a binary system with a mass-ratio biased 
towards unity, the probability that the mass of the most massive star in the cluster changes 
due to the ejection of the initially most massive star can be as large as 20 per cent. Stellar 
collisions increase the maximum-stellar-mass in a large number of clusters when clusters are 
relatively dense (Mod ^ 10'^ Mq and ro.5 = 0.3 pc) and binary-rich. Overall, we conclude 
that dynamical effects hardly influence the observational maximum-stellar-mass - cluster- 
mass relation. 
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1 INTRODUCTION 

Weidner, Kroupa & Bonnell (2010) compiled from the literature 
observational data of 100 young star clusters, whose masses lie be- 
tween ~ 10 and ~ 2 X 10^ Mq and whose ages are younger than 
4 Myr. They showed that observed young star clusters exhibit a 
well-defined correlation between the maximum stellar mass in the 
cluster, mmax, and the mass in stars, i\fcci, of the cluster. An upper 
age limit of 4 Myr was chosen in order to minimize any evolu- 
tionary effects on the sample. The examples of evolutionary effects 
discussed in their paper are as follows. Firstly, mass loss of massive 
stars due to stellar evolution may influence the cluster mass. Sec- 
ondly, gas expulsion leads the cluster to lose a significant amount 
of its stars (i.e. cluster mass) by weakening of the gravitational po- 
tential when the residual gas is expelled from the cluster. However 
these effects unlikel y affect Meci o wing to the young ages of the 
clusters in the Weid ner et al. I( l2010l) sample. The authors corrected 
T7^max for Stellar evolution in the case of O-type stars (note that later 
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than O-type stars would not have evolved much at this young age), 
so that the mmax values provided in their paper can be considered 
as initial values. One process they did not take account of is the 
dynamical ejection of the most massive star from the cluster. The 
authors commented that it is highly unlikely to happen. But this has 
not been studied thoroughly so far. Thus it is our aim in this study 
to investigate how often a young star cluster ejects its most massive 
member. 

This observ ed correlation fits a semi-analytica l model well 
dWeidner & Krou pa 2004, 2006:I Weidner et al.ll2010l) which is de- 
duced from there being exactly one most massive star in the cluster, 



(1) 



^(m)dm, 
subject to the normalisation 

Meci = / m^(m)dm, 



(2) 



where mmax* ~ 1 50 — 300 M,ti is the fundament a l upper limit 
of stellar masses (jWeidner & Kroupal |2004 bOOd iFigetl 120051: 
lOev & Clarkel l2005: 'Crow ther et allboiOl) , min,„ = 0.08 Mq is 
the hydrogen burning mass limit (brown dwarfs contribute neg- 
ligibly to the cluster mass, Thies & Kroupa 2007) and ^(m) is 
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the stellar initial mass function (IMF). A pure size-of-sample ef- 
fect as expected from random sampling has been excluded as 
an origin of the observed correlation. Details of previous studies 
on the mmax-A/cci relation can be found in IWeidner & Kroupal 
( |2004 [200 6). Weidner et al. (2010) and references therein. Note 
that lMaschberger & Clarke! ( 1200 8) argued that at least for the low- 
clusters (i.e. low mass clusters) observed maximum stellar 
masses do not much deviate from random drawing. This is basi- 
cally true but leads the reader to the misinterpretation that a phys- 
ical ^rigin_oftii£inost-massive star in these clusters is ruled out. 
But[Weidner et al.' ('2010) show that in the mass regime considered 
bv lMaschberger & Clarke (2008.) a physical mmax and a stochastic 
mmax ca n not be distinguished from each other. The observed clus- 
ters from IWeidner et al.l (1201 Ol) and their semi-analytical relation 
are reproduced in Fig. 1. Interestingly, a similar relation appears in 
numerical simulations of star cluster formation as well. Smoothed 
particle hydrodynamics numerical simulations of massive star for- 
mation driven by competitive accretion (Bonnell, Vine & Bate 
2004) showed that the most massive star in the forming cluster 
grows following the relation mniax(t) oc Mcci(t)^^^ with time, 
t, which is a best-fit to their simulation data. This fit agrees with 
the semi-analytical mmax-Mcd relation very well for clusters with 
Afeci < 10^ ^0- Furthermore, fragmentation-induced starva- 
tion studied with radiation-hydrodynamical simulation s of massive 
star f ormation using the adaptive-mesh code FLASH I Peters et aEI 
l20IOh also reproduce the relation found by iBonnell et al.l ( 1200411 
Data from the numerical simulations of star cluster formation in- 
cluding the two studies mentioned above are also plotted in Fig. 1. 
The simulation data are in good agreement with the observed data. 

Although most of the clusters follow the relation well, 
there is a spread of mmax values at a given cluster mass. Is 
this spread due to stochastic effects that occur during the for- 
mation of a cluster, or does it mask a true physical functional 
dependence of mmax on Mcci? Dynamical processes can exert 
an influence on the relation during the early evolution of the 
cluster. Stars can be dynamically ejected through energetic few 
body interactions. The lightest star among the interacting stars 
generally obtains the highest velocity and it is unlikely that 
the most massive star is ejected. Several theoretical studies, 
nevertheless, have shown that massive stars can be dynamically 
ejected under certain circumstances such as fro m a small group 
of m assive stars lacking low-mass stars d Clarke & Pringlel 
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through bina ry-single (Hills & FuUertpnl 1980l) and binarv-binarv 
interactions dLeonard & Duncan 1990). The high efficiency of dy- 
namical ejections from dense stellar systems can indeed explain the 
difference between the observed and expected number of OB-type 
stars in the Orion Nebula Cluster dPfl amm-Altenburg & Kroupal 
I2OO6I) . Furthermore, dense and massive RI36-type clusters are 
efficient in expelling massive stars (Banerjee, Kroupa & Oh 2012). 
Thus it may be possible that the heaviest star in a cluster be 
dynamically ejected from the cluster. 

In this contribution we assume there exists an exact function, 
mmax = fn(Mcci), and we study the ejection of the heaviest star 
in a cluster using direct N-body integration to investigate the effect 
on the mmax-Meci relation. Details of the initial conditions of the 
cluster models and of the calculations are described in Section 2 
and then results are shown in Section 3. The discussion and the 
conclusions follow in Sections 4 and 5. 
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Figure 1. Ma ss of the most massiv e star versus cluster mass from observa- 
tional data in lWeidner et al.l l l2010h . Each grey filled circle is a star cluster 
from their table B 1 . Black solid line is the semi-analytical rvimax-Afccl re- 
lation fro m equations 1 and 2 assuming an upper limit of the stellar mass of 
150 Mq iWeidner & Kroupal2004ll2006l) . The open circles are the mass of 
the most massive member in the group vs the total group mass in 14 young 
stellar groups in Taurus, Lupus3, Chal, and IC348 ( Kirk & MverslbOlTh . 
These young (low-mass) stellar groups also follow the rrimax-Afoci rela- 
tion well. Black points and errorbars represent the average and standard de- 
viation of logj^Q Moci and log^o "^max in each bin each of which contains 
22 clusters for the upper 5 bins and 4 clusters for the lowest cluster mass 
bin. Data from several numerical simulations of stai' formation (Bonnell, 
Vine & Bate 2004; Bate 2009, 2012; Smith, Longmore & Bonnell 2009; 
Peters et al. 2010) are included in the figure as open triangles. 



2 MODELS 

We perform a large set of direct N-bod y calculations of young star 
clusters using NBODy6 ( lAarsethlll999l) with various initial condi- 
tions. Cluster masses range from 10 Mq to 10'^'^ Mq with an 
interval of 0.5 on the logarithmic scale and each mass is initialized 
with two different half-mass radii, ro.s = 0.3 pc and 0.8 pc. To 
study the effect of binaries, we adopt two extreme binary fractions 
which are (all stars are single) and I (all stars are in binary sys- 
tems). The initial binary population used in this study is described 
in Section 2.1. Single star clusters (all stars are single) are chosen 
for comparison purpose only since most s tars in actuality form in a 
binary system j Goodwin & KrouDall2005b . 

For each cluster mass the number of stars, A^star, is assigned 
by dividing the cluster mass by the average stellar mass of the clus- 
ter 



iVstar = 



Mo. 



< m > 

where the average stellar mass of the cluster, < m >, is 



(3) 



< m > = 



We adopt the cano nical two-part power law IMF dKroupallioOll; 
iKroupa et alj2012l) . 
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^(m) cx m 
where 



(4) 



a.2 



1.3, 
2.3, 



0.08 s; m/M0 < 0.50, 
0.50 ^ m/M0 ^ 



We use mmax.wK = n^max, which is calculated using the semi- 
analytical relation (equations 1 and 2) assuming the fundamental 
upper limit of stellar masses to be mmax* = 150 Mq. The solid 
line in Fig. 1 represents mmax.wK. 

Individual masses of all stars but one (i.e. A'^star — 1 stars) 
in each cluster are randomly drawn from the IMF (equation 4) 
with a stellar mass range from 0.08 Mq to mmax.wK. To sim- 
plify the analysis, one mmax.wK star is added so that every cluster 
has at least one star with a mass of mmax, wk- El This procedure 
removes the stochastic effects on the initial rrimax. This choice, 
however, gives a bump at the most massive mass-bin in the IMF 
of the cluster, which is especially significant in small- A/^ (i.e. low- 
mass) clusters. As the dynamical ejection of massive stars occurs 
by close encounters betwe en massive stars dClarke & Pringlell992l : 
[Leonard & Duncanlll990h . our enforcement of having a mniax,WK 
star may enhance the ejections at a given Med by overpopulating 
massive stars. We find that this choice would not change our con- 
clusion (see further discussion in Section 4). 

Positions and velocities of stars in the cluster are generated 
according to the Plummer model ( lAarseth. Henon & Wielenll974l) 
which is the si mplest stationary solu tion of the Collisionless Boltz- 
mann Equation (Heggie & Hut'2003 UKroupal2008l ) an d is an excel- 
lent d escription of the nearest star cluster, the Hyades dRoser et al.l 
l20Ilh . To study the effect of initial mass segregation, half of our 
cluster models are initially mass segregated. The method for con- 
structing initially mass segregated clusters in which positions and 
velocities are dependent on stellar masses is described in Section 
2.2. For unsegregated clusters positions and velocities are assigned 
to stars independently of their masses. 

Dynamical time scales such as the crossing time and the me- 
dian two-body relaxation time are important tools to estimate the 
dynamical evolution of stellar systems. The initial crossing time is 

= (5) 



where a is initial velocity dispersion, a — \/ GM^ci Ir 



I gray ~ 2.6 Tp.s i s the gravitational radius dSinnev & Tremain3 
ll987l : lKrouD3l20iD8t) . The relaxation time is 



trel =0.1' 



(6) 



" In iVstar 

Initial conditions of all cluster models are listed in Table 1. 
Table 2 shows the physical properties of the 6 different-mass clus- 
ters. All clusters are evolved up to 5 Myr. And stellar evolution is 
taken into account using the stellar evolution library (Hurley, Pols 
&, Tout 2000) in the nbody6 code. We carry out 100 computations 
for each set of initial conditions. In total 7200 models are thus cal- 
culated with a standard PC. In addition, we perform 10 calcula- 
tions for clusters with A/eci = 10* Mq and with the same initial 
conditions as MS30P in Table 1. We compute only the most ener- 
getic cluster model (MS30P) in our library for Af^ci = lO" Mq 
because of the high computational cost for massive clusters. This 

^ This pr ocedure is not requ ired if optimal sampling of stellar masses from 
the IMF i Kroupa etai] |2012) were used instead of random sampling. Op- 
timal sampling was not available though at the time the present library of 
clusters was computed. 



Table 1. Initial conditions of cluster models. Model name is in the first 
column. Column 2 indicates the initial mass segregation, N standing for 
an initially unsegregated cluster while Y signifies an initially segregated 
cluster Columns 3 and 4 present the initial half-mass radius of the cluster, 
ro.5, and the initial binary fraction, /bi^ respectively. The OP and RP 
in the /bin.i column represent the pairing method for the massive binaries: 
ordered pairing, and random pairing, respectively. The description of the 
pairing methods can be found in Section 2.1. 



Model 


mass segregation 


ni.s [pc] 


/bin,i 


NMS3S 


N 


0.3 





NMS3RP 


N 


0.3 


1 (RP) 


NMS30P 


N 


0.3 


1 (OP) 


NMS8S 


N 


0.8 





NMS8RP 


N 


0.8 


1 (RP) 


NMS80P 


N 


0.8 


1 (OP) 


MS3S 


Y 


0.3 





MS3RP 


Y 


0.3 


1(RP) 


MS30P 


Y 


0.3 


1(0P) 


MS8S 


Y 


0.8 





MS8RP 


Y 


0.8 


1 (RP) 


MS80P 


Y 


0.8 


1 (OP) 



is the largest currently existing systematically generated library of 
young star cluster models. 



2.1 Primordial binaries 

To set up the primordial binaries we require their initial orbital pa- 
rameters such as periods, eccentricities, and mass-ratios. T he initial 
period distribution adopted in this study is equation (8) in lKroup3 
(1995b), 



fp = 2.5 



logio P - 1 



45 + (logio P - ly 



(7) 



where period, P, is in days. With this distribution function, min- 
imum and maximum log-periods, log^Q Pmin and logj,, Pmax, 
are I and 8.43, respectively. This function shows a flat distribution 
at long-periods and is in good agreement with the period distribu- 
tion of low-density young stella r aggregates such as Taurus - Auriga 
iKroupa & Petr-Gotzens 201 ll: iMarks" Kroupa & Od 12011 ). The 
initial eccentr icity distributi on follows the thermal distribution, 

i 'je) = 2e dKroupal l200i) . Pre-main-sequence eigenevolution 
KrouDjl995bl) is not included in our calculations as our emphasis 
is on the massive stars. 

Initial mass-ratios of low-mass stars such as G-, K- 
and M- dwarf b i naries can b e well described with random 
pairmg dKroupal Il995al . |2008|) . Observational studies of OB- 
type binaries, however, show that they tend to have similar 
mass comparuons lOa rcia & Meniiilliod 2001; Sana et al. 200^; 
ISana. Gosset & Evansrr2009l) . ISana & EvansI ( 1201 ll) show that the 
mass-ratio distribution of O star binaries seems uniform in the 
range 0.2 ^ m^lmx ^ 1.0, m.\ being the primary and the 
secondary mass. In any case the mass-ratios of massive binaries 
are high compared to low-mass stars. Random pairing cannot pro- 
duce the observed mass-ratio distribution of massive binaries since 
it typically leads to massive stars being paired with low mass stars 
which are the majority in the cluster. Thus random pairing over the 
whole stellar mass-range for OB star primaries is ruled out by the 
observation. A different pairing method for masses of binary com- 
ponents is needed to create massive binaries. 
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Table 2. Characteristics of clusters with different sizes and masses. Cluster mass, M^^i, number of stars in a cluster, A^star, initial mass of the most 
massive star in the cluster, mmax.i (from equations 1 and 2), and the tidal radius, r^^, are presented in columns 1-4. The rtid is obtained from 
''tid = Rcci^ecl/ (3Mgai))l/'' by assuming Afg^i, the Galactic enclosed mass within the galactocentric distance (i?Gc) "f 8.5 kpc, to be 5 X 10^'' Mq. 
The initial crossing, tcr, and the initial relaxation time, t^ei, for two different cluster sizes, ro.s = 0.3 and 0.8 pc, are given in columns 5-8. Each cluster 
model in Table 1 contains all six different mass clusters in this table. In total thus 72 different initial cluster configurations are used in this study, whereby 100 
random realisations of each configuration are computed with Aarseth's NBODY6. The last model (Afcd = 10^ ^^0) computed 10 times and for the 

MS30P configuration. 



Mecl [Mq] 




mmax,i [Mq] 


»"tid [pc] 


tcr[Myr] 




^rcl 


[Myr] 










{ro.5 = 0.3pc 


O.Spc 


0.3pc 


O.Spc) 




28 


2.1 


3.5 


2.56 


11.13 


2.14 


9.35 


101-5 


76 


4.5 


5.1 


1.44 


6.26 


2.52 


10.98 


102.0 


214 


9.7 


7.4 


0.81 


3.52 


3.22 


14.03 


102-5 


618 


21.2 


10.9 


0.45 


1.98 


4.37 


19.03 


103.0 


1836 


43.9 


16.0 


0.26 


1.11 


6.24 


27.19 


103.5 


5584 


79.2 


23.5 


0.14 


0.63 


9.30 


40.50 


1040 


17298 


114.7 


34.5 


0.08 




14.32 





In this study we introduce a simple method to generate mas- 
sive binaries having mass-ratios biased towards unity. First, all stel- 
lar masses are randomly drawn from the IMF, and then, stars more 
massive than 5 Mq are sorted with decreasing mass and the oth- 
ers are retained in random order. We pair stellar masses in order so 
that a massive star has the next massive one as a companion, while 
stars less massive than 5 Mq have a companion which is randomly 
distributed. Thus binaries with primary masses more massive than 
5 Mq have mass-ratios biased towards unity. We call this method 
"ordered pairing" (OP). Note that star clusters with masses < 100 
Mq contain no stars more massive than 5 Mq in this paper, thus 
OP clusters with M^ci ^ 10^ '^ Mq are the same as random pair- 
ing (RP) clu sters. For a deep discussio n on pairing methods for 
binaries, see iKouwenhoven et al] ^20091) and Weidner, Kroupa & 
Maschberger (2009) who study several pairing mechanisms. 



clusters we use the method introduced in Baumgardt, de Marchi 
& Kroupa (2008). Details of setting up the segregated cluster are 
described in their Appendix. With this method, the heaviest star is 
most bound to the cluster and is located in the core of the clus- 
ter. And the cluster is initially in virial equilibrium and follows 
the Plummer density profile. Although one can vary the degree of 
mass segregation with this method, for simplicity, mass-segregated 
model clusters in this study are fully segregated. Thus the segre- 
gated clusters have a core of massive stars in the centre of the clus- 
ter at the beginning of cluster evolution. 

We stress that such N-body models of initially fully mass- 
segregated clusters with a 100 per cent binary fraction and a mass- 
ratio near unity for the massive binaries have never been attempted 
before. 



2.2 Primordial mass segregation 

Many young star clusters exhibit evidence for mass segregation 
jGouliermis et al.ll2004lchen. de Griis & Zhaoll2007h . It has been 
under debate whether the observed mass segregation of young star 
clusters is the outcome of the star formation processes or of dy- 
namical evolution of the clusters since a certain time is needed for 
it to occur dynamically. Some observed clusters seem too young 
for mass-segregation to have occuiTed. The dynamical mass segre- 
gation time scale is 

tms ~ trel, (8) 

^^massivc 

where mmassivc is the mass of the massive star. For some clusters, 
tms could be shorter than or comparable to their age (e.g. the Orion 
Nebula Cluster which has tms of about 0.1 Myr, Kroupa 2002). 
Thus it is difficult to determine whether an observed mass segre- 
gation is primordial or the result of dynamical evolution. To study 
the influence of primordial mass segregation on the early dynam- 
ical evolution of clusters would give a hint for an answer to this 
problem. However it is beyond the scope of this paper as a deeper 
study on individual clusters is required to do that. 

It is expected that initially mass segregated clusters ought to be 
more efficient in ejecting massive stars thus allowing the distribu- 
tion of massive stars to be used as a constraint on the issue of initial 
mass segregation ( [ Clarke & P ringle 1992; Gvaramad ze & Soman si 
l2008l ; lGvaramadze et al.li201lh . In order to create mass-segregated 



2.3 The A^-body code 

NB0DY6 is a fully collisional A'^-body code which calculates the 
force on a particle from other particles with direct summation. It 
uses the Hermite scheme for integrating the orbits of stars. The 
code adapts individual time steps depending on the local environ- 
ment and the Ahmad-Cohen neighbour scheme jAhmad & CohenI 
1 19731) for calculation efficiency. For treating close encounters, 
Kustaanheimo-Stiefel two-body regularization and Chain regular- 
ization for higher order multiple systems are used. 

Stellar evolution is implemented in the code using fitting func- 
tions with the Single Star Evolution package ( Hurlev et al. 200^) 
and the Binary Star Evolution package (Hurley, Tout & Pols 2002), 
which allows a collision between two components of a binary. De - 
tails of stellar evolution in nbody6 can be found in lHurlevl ( l2008h . 
Since we activate stellar evolution in the code, a star has a radius 
instead of being a point mass particle and may collide with other 
single stars or its companion in a binary system by close encoun- 
ters and/or binary hardening. When they merge, a mass of one star 
is replaced by the sum of the colliding stars and the other star is 
replaced by a massless particle with a large distance so that it is 
removed from the calculation as a massless escaper. A metallicity 
of Z = 0.02 (solar) is adopted in all our calculations. 
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3 RESULTS 

Clusters keep losing their mass due to stars escaping from them via 
the two-body relaxation process or through dynamical ejections be- 
sides stellar evolution, and so we only count stars found inside the 
tidal radius as cluster members. A cluster mass does not change 
much (at most « 3 per cent on average by 3 Myr) within the first 
few Myr. Therefore the change of the mmax-Mcd relation within 
the first few Myr is mainly caused by the change of mmax- There 
are three ways to change the maximum stellar mass as a cluster 
evolves. First, stellar evolution changes the mass of the heaviest 
star in the cluster. Stars lose their mass with time via stellar winds. 
The mass loss rate is dependent on the stellar mass. The more mas- 
sive a star is the larger its mass-loss rate is. In the case of the most 
massive cluster in our model (Afcd = 10"^'^ Mq), its heaviest star 
with an initial mass of ~ 80 Mq loses mass to become a « 64 Mq 
star at 3 Myr On the other hand for clusters with i\feci ^ 10^ Mq, 
Timax remains almost the same as the initial value since the heav- 
iest stars in such clusters do not evolve much in a few Myr. Sec- 
ondly, stars can physically collide. If two stars collide and become 
a more massive star than the initially heaviest star, or if the initially 
heaviest star collides with another star, then the new heaviest star 
will lie off the initial relation. Massive stars generally move into the 
cluster centre which has a high stellar density. Thus a massive star 
may collide with another massive star in the cluster centre. Lastly, 
dynamical ejection of the initially heaviest star in the cluster also 
changes mmax by replacing it with the initially second massive star 
in the cluster. 

Stellar evolution and dynamical ejection leads to the cluster 
having a smaller maximum stellar mass, while stellar collisions in- 
crease the maximum stellar mass in the cluster. In our models all 
three of these effects occur. However here we concentrate only on 
stellar collisions and ejections as the initial masses of the heaviest 
stars are set to be equal for the same cluster mass. Thus stellar evo- 
lution does not produce the mmax spread at the same cluster mass. 
And, we are particularly interested how often the heaviest star is 
ejected from the clusters. 

Using direct N-body calculations, we study the dynamical ef- 
fect on the mmax-Mcci relation with various initial conditions of 
the clusters (Tables 1 & 2). Although the clusters are evolved up to 
5 Myr, we only use the results up to 3 Myr since stellar evolution 
begins to play a dominant role in changing the value of mmax of 
the lO'^'^ Mq cluster at around 3.5 Myr. The initially most mas- 
sive star of the 10^'^ Mq cluster becomes a blackhole after 4 Myr 
Furthermore, most of the observed cluster s used in the study of the 
mmax-Afcci relation dWeidner et al.l2oiol) are younger than 3 Myr. 

For clarification, we refer to the heaviest star in the cluster at 
Myr as Smaxi and to the mass of the most massive star in the 
cluster at a given snapshot as mmax- 

3.1 Dynamical ejection of Smaxi 

Stars can be ejected from a cluster via close encounters between 
a hard binary and a single star/binary. During the encounter, the 
hard binary gives its binding energy to the interacting star/stars 
in the form of kinetic energy and it hardens. The star that gained 
kinetic energy may be ejected from the cluster. Generally the 
lightest one among the interacting stars attains the highest veloc- 
ity after the interaction. Thus the dynamical ejection of massive 
stars preferably occurs from interactions between massive stars. 
[Leonard & Dunca n ( J_990) showed that binary-binary interactions 
are the most efficient way for producing massive runaways. 



We consider Smaxi to be dynamically ejected if the star is 
further away from the cluster centre than the tidal radius of the 
cluster. The number of clusters which eject their Smaxi, A'rcsc, is 
listed in Table 3 for all cluster models. 

In the following subsections, results on the dynamical ejection 
of Smaxi are discussed separately for models with different initial 
half-mass radii. 



3.1.1 The ro.5 = 0.8 pc models 

In Table 3 there are two clusters with ro.s = 0.8 pc whose Smaxi is 
located further than the tidal radius of the cluster at 3 Myr. The most 
massive star of one not initially mass-segregated 10 Mq cluster 
was, in fact, located outside of the tidal radius of the cluster, which 
is ~ 3.5 pc, at Myr. Therefore this case is not due to dynamical 
ejection. Thus only I out of 3600 model clusters with ro.s = 0.8 pc 
eject their Smaxi within 3 Myr. 

Fig. 2 shows the mmax-A/cci relation of the binary-rich clus- 
ters with ro.5 = 0.8 pc and massive binaries paired by the OP 
method (models NMS80P and MS SOP) at different ages of the 
clusters. Even though this set of initial conditions is the most dy- 
namic case of the cluster models with ro.s = 0.8 pc, only one 
cluster, with a mass of 10^'^ ejects its Smaxi (Table 3). But 
the initially second heaviest star of the cluster, which becomes the 
most massive one in the cluster after ejection of Smaxi, has a sim- 
ilar mass to the mass of Smaxi. Therefore the effect of the ejection 
on the mmax-.Aioci relation is negligible in this case. 

It is unlikely that the mmax-i^ioci relation is affected by dy- 
namical ejection of Smaxi for clusters of this size. But it is worthy 
to note that at 3 Myr a few OP clusters show their Smaxi moving 
faster than the escape velocity, Vcadj') = \/2|$(r) |, where <l>(r) 
is the gravitational potential at a distance r from the cluster centre, 
although these clusters barely eject their Smaxi (Table 3). 

3.1.2 The ro.5 = 0.3 pc models 

Fig. 3 shows the mmax-Mcci relation for the binary-rich cluster 
models with ro.s — 0.3 pc. Despite the smaller size of the clus- 
ters, none of the single star clusters (NMS3S and MS3S) eject 
their most massive star (Table 3). Only two massive clusters with 
Med ~ 10'^'^Mq eject their Smaxi when the massive binaries 
are randomly paired (NMS3RP and MS3RP in Table 3). However 
binaries help the most massive star to attain a higher velocity com- 
pared to the single-star clusters. The number of clusters for which 
the heaviest star has a speed exceeding the escape velocity is larger 
when the stars are initially in a binary system (Table 3). 

Clusters with Afcd ^ 10^'^ Mq hardly eject their most mas- 
sive star even though the clusters form in energetic initial con- 
ditions such as being binary-rich and mass-segregated. Among 
all clusters with Med ^ 10^'^ Mq, only one out of 100 clus- 
ters with Mcci = 10^ Mq and two out of 100 clusters with 
Meci = lO^ '^ M0 eject their Smaxi (Table 3). However at 3 Myr 
there are a few clusters whose Smaxi is inside the tidal radius but 
has a velocity higher than the escape velocity so that it may leave 
the cluster later. 

For clusters with il/ed ^ 10^ -^0' 15-40 per cent of the most 
energetic models (MS30P) have Smaxi moving faster than the es- 
cape velocity and 2-20 per cent of the clusters have ejected their 
Smaxi at 3 Myr. When the clusters are initially mass-segregated, 
slightly more clusters eject Smaxi. 

Fig. 4 shows the ejection frequency of Smaxi as a function 
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Figure 2. The mass of the most-massive star versus the cluster mass for the unsegregated (left: NMS80P) and the segregated (right: MS80P) clusters with 
''0.5 = 0.8 pc and massive binaries with mass-ratios biased towards unity at 0, 1, 2, and 3 Myr (from top to bottom). Each black point indicates a cluster The 
grey data are the average observational data as in Fig. 1. The solid line represents the semi-analytical model from Weidner & Kroupa (2004, 2006) assuming 
an upper hmit of stellar mass of 150 Mq (equations 1 and 2). As we fix the number of stars at a certain cluster mass, the cluster mass slightly varies for each 
realization for the same cluster mass model. We only plot the result up to 3 Myr since the initially heaviest star of the cluster with 10^'^ Mq loses a large 
amount of its initial mass within 3.5 Myr, thereafter stellar evolution affects the relation. The stars that appear above the solid curve are merger products. 



of cluster mass for the MS30P models. We only plot these models 
since the other models barely eject their initially most massive star 
(see Table 3 for all models). The ejection probability of the Smaxi 
star increases with the cluster mass as the stellar density increases. 
About 8 (20) per cent of 10^ '^ (10*) Mq clusters eject their Smaxi 
within 3 Myr (squares in Fig. 4). Some massive clusters which eject 
their Smaxi can be missed as we use the tidal radius as a criterion 
for the ejection. Some ejected Smaxi from massive clusters may 
not reach the clusters' tidal radius by 3 Myr due to their large tidal 
radius. Thus, the real ejection frequency would be higher than the 
above value for the massive clusters. By using A'^vcscio (Table 3), 
the number of clusters whose Smaxi has a velocity greater than the 
escape velocity of the cluster and has travelled beyond 10 pc from 
the cluster centre, the ejection frequency increases to 40 per cent for 
the cluster with 10* Mq (circles in Fig. 4). In order to estimate the 
real probability for the Smaxi ejection, in addition, we provide in 



Appendix A the number of clusters whose Smaxi is found beyond 
the distance criterion using the half-mass radius for all models. 

From Figs 2 and 3 it is evident that collisions of massive stars 
may conceal the dynamical ejection of Smaxi by the product of 
the collisions becoming more massive than the mass of the initially 
heaviest member. Thus the dynamical behaviour of the Smaxi 
needs to be considered to distinguish whether it is ejected or not. 
Figs 5 and 6 show the distance from the cluster centre and the ve- 
locity of the Smaxi in NMS30P and MS30P clusters with 10 Mq 
and 10^'^ Mq. 

The heaviest stars are initially located at a wide range of radii 
up to ~ 2.5 pc in the unsegregated cluster models (Fig. 5) while 
they are centrally concentrated in the segregated ones (Fig. 6). In 
both cases, massive stars sink towards the centre of the clusters 
due to dynamical friction and/or energy equipartition. This is more 
prominent in the unsegregated clusters since the heaviest stars al- 
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Figure 3. Same as Figs 1 and 2, but for the OP clusters with ro.s = 0.3 pc {left: NMS30P, right: MS30P). Note how including OP binaries and mass 
segregation increases the scatter in the diagrams. 



ready reside in the deep potential of the segregated clusters at the 
beginning of the calculations. 

As shown in Figs 5 and 6, the heaviest stars hardly attain a 
high velocity if the massive stars are randomly paired into binaries. 
In random pairing, massive stars are mostly paired with low mass 
stars therefore their mass ratios (m2/mi, where mi m2) are 
skewed to 0. Randomly paired massive binaries therefore behave 
like single stars. Clusters with massive binaries paired randomly do 
not shoot out their Smaxi more frequently even though the clusters 
are initially mass-segregated. 

Clusters with massive binaries paired by OP and Med ^ 
100 Mq effectively produce heaviest stars with velocities exceed- 
ing the escape velocity of the cluster For example, in the case of 
some initial conditions more than 20 per cent of the clusters show 
that their Smaxi has a velocity larger than the escape velocity at 
3 Myr (Table 3). It is known that massive stars are ejected from 
the small core of massive stars in the cluster centre lacking low 
mass stars ( IClarke&Pringlelll992h . In the case of segregated clus- 
ters with OP massive binaries, the clusters already form with this 
kind of core thus the massive stars can be ejected at a very early 
age of the cluster. S'maxi of one cluster with A/eci ~ 10'^'^ Mq 



from the MS30P model has travelled more than 200 pc from the 
cluster with a velocity of about 80 km at 3 Myr (Fig. 6). 

Although the low mass clusters barely eject their Smaxi re- 
gardless of their initial conditions, the binary-rich or binary-poor 
cases show differences, e.g. the number of clusters whose Smaxi 
has a velocity greater than the escape velocity is larger for binary- 
rich clusters. For massive clusters, on the contrary, binary-poor and 
binary-rich clusters with random pairing show similar results in 
the case of unsegregated models. Primordial mass segregation en- 
hances the ejection of the most massive star at earlier times and 
helps more clusters shoot out their Smaxi within 3 Myr 

3.2 Stellar collisions 

mmax can increase by stellar collisions, either involving Smaxi 
or not, which makes a star heavier than the mass of Smaxi be- 
fore the collisions occur Although stellar collisions occur over a 
whole range of stellar masses, in this study, we only care about the 
collision that changes mmax. The occurrences of the collisions are 
contained in Table 3. 

For clusters with ro.s — 0.8 pc (2 Mq pc"'^ < po.5 < 
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Figure 5. Distances from the cluster centre and velocities of Smaxi stars of the initially unsegregated binary-rich clusters with ro.5 = 0.3 pc {left: NMS3RP, 
right: NMS30P). Each dot denotes a cluster and in total there are 100 dots per configuration. Squares and circles are the values of Smaxi in the clusters with 
10 Mq and 10^'^ -^0> respectively. Filled (red) symbols represent Smaxi stars with a speed exceeding the escape velocity. The black (grey) solid curves 
are the escape velocity of one cluster with 10'^'^ (10) Mq as a function of distance from the cluster centre at each Myr The grey dotted and the black dashed 
vertical lines indicate the initial tidal radius of the cluster with 10 Mq and with lO'^'^ '^O- respectively. Note how mass segregation develops by 1 Myr and 
how OP increases the occurrence of ejected most massive stars. 



750 Mq pc~'^, where po.5 = 3A'/oci/(87rro.5) is the average mass 
density within the half-mass radius) stellar collisions rarely occur in 
them due to their low density. In single-star cluster models (NMS8S 
and MS8S), only one cluster shows that its mmax changes by a 
stellar collision which involves Smaxi- In the case of binary-rich 
models there are a few clusters whose mmax changes via stellar 
collisions, mostly in the most massive cluster models, but with a 
probability of less than 10 per cent taking all models into account. 
In low-density (ro.5 = 0.8 pc) clusters a change of mmax through 
stellar collisions is highly improbable. 

For clusters with ro.5 = 0.3 pc, the stellar collision result of 
low-mass clusters, A/od lO^ '^ Mq (po.5 < 1400 Mq pc"^) 
is similar to the result from clusters with ro.5 = 0.8 pc. But in the 
case of the massive cluster models (A/cd 10"^ Mq), especially 
the most massive one, about half of them experience a change of 
mmax by stellar collisions. 

Most of the collisions are induced by binary-encounters, in 
which the two components of a binary system collide due to their 



highly eccentric orbit generated by perturbation through other stars. 
Direct dynamical collisions are extremely rare. This naturally ex- 
plains why TVc in Table 3 becomes larger when clusters are initially 
binary-rich. And the number increases when massive binaries are 
paired with the OP method when compared to RP. This can be un- 
derstood because both components of a massive binary in the OP 
method are massive stars and thus have larger sizes leading to col- 
lisions. With a different collision channel from our result, Gaburov, 
Gualandris & Portegies Zwart (2008) also showed that a binary is 
more efficient in stellar collisions than a single star as the enhanced 
cross-section of a binary compared to a single star results in other 
stars engaging the binary, and then this can lead to a collision be- 
tween one of the binary components and the incoming third star 

Both the dynamical ejection of the Smaxi and stellar col- 
lisions that increase mmax barely take place in the same clus- 
ter w ithin the cluste r mas s ra nge we study. As show n by 
iBaumga rdt & Klessert j201lh and lMoeckel & Clarkd jlOllI) . stel- 
lar collisions generally lead to the formation of a single very mas- 
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Figure 6. Same as Fig. 5 but for the initially mass-segregated binary-rich clusters with rp.s = 0.3 pc (left: MS3RP, right: MS30P). 



sive star through the merging of several massive stars rather than 
the formation of many massive stars. This single very massive star 
is hardly ejected from the cluster and its formation reduces the 
number of massive stars thus it may hamper the ejection of other 
massive stars. Out of 7200 cluster models (excluding the 10* Mq 
clusters from MS30P model), there are only three clusters in which 
both events occur. In one cluster Smaxi is dynamically ejected, but 
one massive binary in the cluster merges then becoming more mas- 
sive than Smaxi- In the other two clusters their Smaxi gains mass 
by the collision with another star and then it is dynamically ejected. 

Although this event is very rare, the most massive star 
in the star-forming region LH 95 in Large Magellanic Cloud 
JPa Rio et alj l2012f) might be an example. The peculiarities of 
the star, with a much younger age than the average age of other 
sta rs in the reg ion and the mass being higher than mmax from 
the lWeidner & Kroupa (2004) 7nmax--A^Gci relation, could have re- 
sulted from a stellar collision between binary components with 
similar masses of which the primary star may have been the ini- 
tially most massive star of one of the three main substructures in 
the region, with a mass following the mnmx--A^cci relation from 
IWeidner & Kroupa ([2004). The stellar collision could have in- 
creased the stellar mass and rejuvenated the star. And the ejection 



of the star with a low velocity can explain the location of the star 
which is at a rather far (~ 10 pc) distance from any of the substruc- 
tures but still within the region. 



3.3 The spread of rrimax 

Fig. 7 presents the standard deviation of log mmax, crimmax, for 
the observed samples and our models. The observed larger rrimax 
spread than what emerges from our models may be a result of 
stochastic effects of star formation as the dynamical processes 
hardly influence the change of mmax in these clusters, as shown by 
this study. However, numerical simulations of star cluster forma- 
tion show that mmax and Mod evolve tightly following the relation 
of mmax oc M^'^ dSonnell et aI.ll2004l : |Peters et al.ll2010l) . 

For low-mass clusters (Afcci < 100 Mq) the differences be- 
tween (Timmax of the obsetvatiou and our model are large despite 
taking into account the spread due to the different cluster masses 
in the bins of observational data. This could be due to the large un- 
certai nties of the ob s ervati ons. Note that the homogeneous data set 
from Kirk & Mverj ( I20T ih is more confined to the relation than the 
inhomogeneous data from lWeidner et all ( I2OI0I) which come from 
many different references. 
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Table 3. Results at 3 Myr. M^^i is a cluster mass in Mq (per Afcci value there are 100 clusters, but 10 clusters for 10* Mq). Nresc is the number of clusters 
whose Smaxi is located beyond the tidal radius of the cluster A^vosc , the number of clusters whose Smaxi has a velocity larger than the escape velocity. 
^vesclOi the number of clusters whose Smaxi is located further than 10 pc from the cluster centre and has a velocity larger than the escape velocity. A^c is 
the number of clusters whose mmax changes due to stellar collisions. Numbers in brackets indicate the collision products that do not involve Smaxi- For 
example, in the case of clusters with M^^i = lO'^'^ Mq from the NMS3RP model, 1 out of 100 clusters lost their Smaxi by dynamical ejection, Smaxi of 
6 clusters have a velocity greater than the escape velocity of the cluster, for 2 clusters out of these 6 clusters the star is located beyond 10 pc from the cluster 
centre. Stellar collisions which change the mmax have occurred in 51 clusters, for 1 1 clusters out of these 51 the collisions do not involve Smaxi- 
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The spread in our models for relatively massive clusters 
(Mcci ^ 10"^ ^q) comparable to the observed one in the case 
of clusters with ro.s = 0.3 pc and with OP for massive binaries 
(NMS30P and MS30P models). 

Stellar collisions dominate the change of rrimax of relatively 
dense clusters with Meci ^ 10^'^ Mq and ro.s ~ 0.3 pc (Ta- 
ble 3) and therefore the spread of mmax of these clusters mostly 
comes from the collisions. However it should be noted that the ex- 
act solution for the stellar collision process is as yet poorly un- 
derstood. The treatment of stellar collisions in the code is simply 



adding the masses of two stars. Thus the collision products shown 
in this study provide only a rough idea about merger rates and their 
masses. Therefore we cannot quantify the real spread produced by 
the collisions but we expect that it would be smaller than our result 
since in reality not all the stellar mass ends up being in the merger 
product. 



© RAS, MNRAS 000,[TltT3 



The dynamical influence on the mmax-Msci relation 1 1 



0.8 



0.6 



0.4 - 



0.2 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

MS30P 
t = 3 Myr 



' rej.SMAXI 
vlOej.SMAXI 



■H- 




2 3 4 



Figure 4. Ejection frequency of the Smaxi as a function of cluster mass 
for mass-segregated clusters witli massive binaries (MS30P) at 3 Myr. 
/roj.SMAxi' rnarked witli squares, is tlie ratio of tlie number of clusters 
whose Smaxi is located further than the tidal radius from the cluster centre, 
A^rosc in Table 3, to the total number of the realizations, A^run, which is 100 
in our study for each initial condition set (10 for 10* Mq clusters). EiTor- 
bars indicate Poisson uncertainties. All other models have /rej, Smaxi ~ ^ 
(Table 3). The circles, /viOoj, Smaxi ' °^ ^'^vcscio (Table 3) to 

.'Vrun- /viOoj, Smaxi probably closer to the real ejection frequency as 
some massive clusters which eject their Smaxi can be missed in A^rcsc 
due to their large tidal radii. The grey star is the ejection frequency of 
S maxi from clusters w ith Mgci = 10^ taken from the calculations 
bv lBaneriee et al.l )2012l) . Those authors refer to a star that is ejected when 
its distance from the cluster centre is larger than 10 pc. Note that the initial 
conditions of their calculations are different but compai'able to our MS30P 
models. 



4 DISCUSSION 

Fig. 8 shows the initial- and the final (t = 3 Myr) mass functions of 
MS30P model clusters. Massive stars are overpopulated in a clus- 
ter with a bump at the most massive stellar-mass-bin of the cluster's 
mass function because we enforce all clusters to have a star with a 
mass of mmax.wK- As a necessary condition for massive star ejec- 
tions is a small core of massive stars, the overpopulation of massive 
stars would enhance the ejection of massive stars. The other effect 
that might result from this choice is the over-spread of mmax in 
low-mass (i.e. small-number) clusters if many of the clusters eject 
their Smaxi- When stellar masses are randomly drawn from the 
IMF low-mass clusters would hardly have a star with a mass close 
to rrimax while massive clusters, that populate stars over the whole 
mass range, would have a few. However, the dynamical ejection of 
the most massive star from the low-mass cluster is highly improb- 
able thus enforcing clusters to have a star with a mmax.wK would 
not have an impact on the spread. 

The orbital parameters of the massive binaries are important 
for the ejection of massive stars. But, our knowledge of their initial 
distribution functions is still poor. We use the same period distri- 
bution for massive binaries as that of low-mass stars. However, 
a high fraction of short-period massive binaries is suggested by 
observations ( S ana & E vans 201 1) implying that massive binaries 
may have a different initial period distribution, compared to that of 
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^immax' ^obs' is much larger than 
the dispersion in the models. To quantify the deviation due to binning the 
cluster mass in the observed range of cluster masses in each mass bin, we 
obtain the mmax according to equations 1 and 2 with m'^^y^ = 150 Mq 
for the observed cluster masses. We then calculate the standard deviation 
of logj^Q mmax, o"i, by binning the clusters in the same way as for Cobs- 
This is shown as grey triangles connected with a solid line. The corrected 
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low-mass binaries. Dynamical ejections of the most massive star 
might become more likely than presented here. Future work taking 
into account the most recent constraints on the period distribution 
of massive star-binaries will investigate this issue. The results for 
Mcci ^ 10'' Mq clusters presented here will, however, not be af- 
fected as such clusters do not contain many, if any, massive stars. 

Apart from our MS30P sequence of models, our young star 
cluster library does not contain clusters more massive than 10* Mq 
as they are computationally expensive due to the large binary popu- 
lation, while the observational sample contains a few clusters more 
massive than 10* Mq. At a higher density massive clusters may 
show the ejection of the heaviest star with a higher probability. In- 
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Figure 8. Mass functions of MS30P clusters. Each line indicates the av- 
erage mass function of clusters obtained from the 100 computations per 
models with a mass of 10, lO^ S, 10^, lO^ S, 10^ and lO^ S Mq from the 
bottom to the top. The solid and dashed histograms are mass functions at 
(initial) and 3 Myr. respectively. A solid line on top of the histograms is the 
Salpeter mass function. 



deed, 20 per cent of the MS30P clusters with Afcd = W Mq 
eject their Smaxi (Fig. 4). Furthermore , 3 out of 4 cl usters with 
Afoci = 10^ Mq computed by Baneri ee et al.l ( 1201 2[) eject their 
Smaxi- Further work including such massive clusters will be car- 
ried out to study this issue. 



5 CONCLUSIONS 

We have established a large theoretical young star cluster library 
using Aarseth's direct N-body integration code, nbody6, with 72 
different combinations of initial conditions such as cluster size and 
mass, initial binary population, and primordial mass segregation. 
The library contains two different sizes of clusters, ro.s = 0.3 and 
0.8 pc, with two different binary fractions, zero and unity. In order 
to take into account that observed OB binaries likely have com- 
panions with similar masses, we generate clusters with massive bi- 
naries not only paired randomly from the IMF but also having a 
mass-ratio close to unity. In both cases the stars follow the canoni- 
cal IMF. And we model initially mass-segregated clusters as well as 
unsegregated clusters. Using this library we study the rrimax-Mcci 
relation during the early evolution (^3 Myr) of the star clusters fo- 
cusing on the effects on it established through dynamical ejection 
of the heaviest star from the cluster under the various initial condi- 
tions. Such computations of fully mass-segregated binary-rich clus- 
ters have never been performed before. 

Stellar evolution, stellar collision and dynamical ejection can 
alter the nimax of the cluster. In our models, all three effects are 
observed to affect rriniax- Stellar evolution affects the relation only 
for the massive (> 10^ ^q) clusters since only stars more mas- 
sive than 40 Mq significantly lose their mass within 3 Myr. Fur- 
thermore, since the mass of the Smaxi is the same for the same 



cluster mass, stellar evolution cannot contribute to the spread of 
the rrimax in this study. Stellar collisions influence mainly massive 
clusters, especially for those with the smaller radius as the clusters 
are denser. For clusters with ro.s = 0.3 pc and Mcci = 10"^'^ 
for about half of the clusters their rrimax has changed by stellar col- 
lisions within 3 Myr. Lastly, concerning the focus of this study, the 
dynamical ejection of Smaxi only occurs in the binary-rich clus- 
ters with ro.5 = 0.3 pc and Meci > lO'^ ^Q - The number of the 
ejections increase when massive binaries are paired with the OP 
method and/or are initially concentrated in the cluster centre. As 
massive clusters likely have a few stars with masses close to the 
mass of Smaxi, "imax value would change little for massive clus- 
ters when Smaxi is dynamically ejected. Overall we conclude that 
(dynamical) evolutionary effects hardly produce a spread of mmax 
in the relation for low-mass {M^ci 10^'^ ^q) or 1^^^ dense 
(ro.s = 0.8 pc) clusters. Massive {Mad = 10^'^ Mq), binary-rich 
clusters with ro.s = 0.3 pc do show a significant spread of mmax, 
mostly produced by stellar collisions, becoming comparable to the 
observed spread (Fig. 7). 

Concerning the dynamical ejection of Smaxi, we find that in 
general it is very unlikely even in relatively dense clusters with 
massive binaries with similar mass companions. However, its prob- 
ability can reach up to 20 per cent in our cluster library depending 
on the initial configuration of the cluster. For example, none of the 
clusters without initial mass segregation but with OP binaries and 
with Afoci = lO^ '' Mq and ro.s = 0.8 pc (NMS80P model) eject 
their Smaxi, while 8 (2) out of 100 (10) mass-segregated clusters 
with Meci = lO^ '^ (lO"*) Mq and ro.s = 0.3 pc (MS30P model) 
eject their Smaxi (Table 3). In reality, many young star clusters 
are observed to fulfil the conditions which are needed to eject their 
massive star, such as having a compact size and massive binaries 
with similar component masses, so that some (up to ~ 75 per 
cent for 10^ Mq clusters, see Fig. 4) of the real embedded clus- 
ters could have lost their initially most massive star by dynamical 
ejection within 3 Myr. 

In this paper, we are only interested in the heaviest star 
of the cluster. However the dynamical ejection of other mas- 
sive stars in the cluster is also interesting because it is impor- 
tant to understand the origin of field massive stars and OB run- 
aways jFuiii & Portegies Zwartll2oTll) to help us to constrain the 
initial configuration of massive stars in clusters. The dynamical 
ejections of OB stars in our theoretical young star cluster library 
wil l be discussed in th e fo llowing paper (Oh et al. in prep.) , 
and lSaneriee et alj ( l2012l) and lPflamm- Altenburg & KroupjHooeh 
have, respectively, already demonstrated that R136-type star-burst 
and ONC-type clusters are very efficient in ejecting OB stars. 
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APPENDIX A: EJECTION ESTIMATES OF Smaxi USING 
THE HALF-MASS RADIUS 

The tidal radius is used in the paper to determine if a star belongs 
to the cluster. In some cases, even though Smaxi is dynamically 
ejected through a strong close encounter the star may have not left 
the cluster yet at a given time, especially if it has obtained a rel- 
atively low velocity from the encounter. Thus the real probability 
for ejection of Smaxi may be higher than that presented in this pa- 
per. Considering that the massive stars are generally seated in the 
central part of the cluster due to dynamical interactions, relatively 
small distance such as the half-mass radius can be used as the ejec- 
tion criterion for Smaxi. However, it should be noted that Smaxi 
being found outside of the half-mass radius does not necessarily 
mean that is ejected, especially when the cluster is initially unseg- 
regated and dynamically unevolved. Here, therefore, we provide 
tables (Table Al and A2) which contain the mean half-mass radius 
and the numbers of clusters whose Smaxi is found outside of ro.s, 
2ro.5 and 4ro.5 and additionally has a velocity greater than the es- 
cape velocity at the given time. For most of the cluster models the 
mean half-mass radii do not change much within 3 Myr due to the 
clusters' dynamical evolution. Note that the initial crossing time of 
the clusters with M^d ^ 100 Mq and ro.s = 0.8 pc is longer than 
3 Myr (Table 2). The most expanded half-mass radius is ~ 1 pc at 
3 Myr, although a few massive cluster models have expanded up to 
twice their initial half-mass radius. 

The initial positions of Smaxi distinctly show whether the 
models are generated with initial mass segregation. For the initially 
unsegregated clusters (Table Al), more than half of the realizations 
for each models show that the Smaxi is initially located outside of 
the half-mass radius. Even in up to 23 (8) per cent of the clusters 
the star is located further out than 2ro.5 (4ro.5) at Myr. For the 
initially mass-segregated clusters, only in few low-mass clusters 
(less than 5 per cent at most) is Smaxi located outside of ro.s at 
Myr (Table A2). It may be strange for the mass-segregated cluster 
to show its most massive star having an initial position beyond the 
half-mass radius at all. But for the initially mass-segregated clusters 
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with low-masses (Med < 100 Mq) the initial position of Smaxi 
can, sometimes, be generated slightly beyond the half-mass radius 
due to the shallow gravitational potential, the low number statistics 
and the algorithm generating the initially mass-segregated cluster 
(the heaviest star being the most bound to the cluster, but not re- 
quiring it to be at the most central position). None of the initially 
mass-segregated clusters have Smaxi located initially further out 
than 2ro.5. 

For the initially not mass-segregated clusters with relatively 
high masses 100 Mq) and/or ro.s = 0.3 pc, N{rsMAxi > 
ro.s) in Table Al significantly decrease at 3 Myr as a result of dy- 
namical evolution, i.e. dynamical mass segregation. But there are 
still a large number of the unsegregated clusters, particularly the 
low-density ones, which have Smaxi beyond the half-mass radius. 
In many of the cases the star has a velocity lower than the escape 
velocity. As those clusters are dynamically young and the velocity 
of the star is too small to escape from its cluster, it is unlikely that 
these stars are dynamically ejected. Either the star was initially lo- 
cated beyond the half-mass radius and has not fallen into the cluster 
centre yet. Or, if Smaxi of the low-density clusters has an initial or- 
bit as large as the half-mass radius (note that the half-mass radius is 
less than a pc), it could be temporarily found beyond the half-mass 
radius since the dynamical interactions which lead to the massive 
star being confined to the central part of the cluster are insufficient 
by 3 Myr for these low-density clusters. This can also explain that 
for the mass-segregated clusters with low masses, especially with 
10 Mq, A'^(?'Smaxi > ''0.5) increases at 3 Myr and most of these 
clusters have Smaxi moving slower than the escape velocity. There 
are only two initially mass segregated clusters with 10 Mq whose 
Smaxi is found beyond twice the half-mass radius at 3 Myr. There 
are due to low-energy encounters expelling the stars from the core 
without ejecting them from the cluster. However, in the NMS30P 
and MS30P models, and in massive (^ lO'^ ^q) clusters from 
the other models, a number of clusters have Smaxi being beyond 
the half-mass radius with a velocity larger than the escape velocity, 
which implies that the Smaxi stars are dynamically ejected. The 
probability of Smaxi ejection increases with cluster mass. 40 per 
cent of the 10* Mq MS30P clusters have probably ejected their 
Smaxi. 
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Table Al. The averaged half-mass radius, < ro.s >, and the number of clusters whose Smaxi is found beyond the half-mass radius, ^(^Smaxi ''0.5 )> 
beyond 2ro.5, ^ {"'"S^k-ki ^ 2ro.5), and beyond 4ro.5, ^(''Smaxi ^ ^r'o.s) at and 3 Myr from the initially unsegregated clusters. The latter numbers at 
3 Myr marked with t, are the number of clusters whose Smaxi fulfils the distance criteria and moves faster than the escape velocity ('^Smaxi ''osc). 
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Table A2. Same as Table Al but for the initially mass-segregated cluster models. Note that only 10 realizations are performed for IC* Mq clusters due to the 
expensive computing cost. 
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